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摘 要 : 计算 速度 对 于 期 权 交 易 者 至 关 重 要 ， 关 系 到 如 何 有 效 地 制定 价格 并 评估 相应 的 风险 ， 而 云 并 行 计算 提供 的 随 
收 随 付 制 (pay-as-you-go) 可 以 实现 低 成 本 运行 。 在 微软 云 平台 Windows Azure 的 基础 上 ， 开 发 了 基于 云 并 行 计算 的 
期 权 定价 试点 云 软件 AzureOP， 该 软件 以 较 低 的 费用 提供 了 低 风 险 和 高 速度 ， 并 给 出 了 AzureOP 对 于 美式 期 权 价 格 的 
模拟 结果 ， 绘 制 了 对 应 的 期 权 价 格 定价 曲线 和 定价 曲面 。 最 后 ， 对 云 并 行 计算 在 金融 应 用 上 的 优势 和 不 足 进行 了 总 结 
和 讨论 ， 同 时 举例 说 明了 试点 云 软件 AzureOP 的 具体 细节 。 
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Abstract: The computational speed is important for the crucial for option traders to efficiently decide the price and evaluate the 
corresponding risk, and cloud parallel computing provides a low-cost implementation by providers’ Pay-as-You-Go policy. 
Based on Microsoft’s cloud platform Windows Azure, developed a option pricing’s pilot cloud software - AzureOP, this software 
provided low risk and high speed with a lower cost, and presented simulation results of AzureOP for American Put, plotted the 
corresponding option pricing curve and pricing surface. Finally, makes conclusions and discusses the advantage and the 
disadvantage of computational finance on the cloud, and illustrates the details development of pilot cloud software AzureOP. 
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[0 图 1 所 示 。 此 PaaS 基于 微软 的 数据 中 心 (IaaS )， 逻 辑 结构 
Windows Azure 的 HPC 调度 程序 和 Azure 节点 组 成 。 选 择 
虽然 云 计 算 被 定义 为 计算 和 存储 交付 ， 但 作为 基于 互联 网 好 PaaS 提供 商 后 ，SaaS 得 到 开发 ， 测 试 和 交付 ， 并 且 成 功 发 
的 计算 0 , 云 计算 也 可 以 被 定义 为 金融 公司 和 制造 公司 的 基于  ” 布 的 SaaS 包括 生物 信息 学 [49， 图 像 处 理 079， 工 程 34， 医 学 信 
计算 的 服务 器 : 计算 密集 的 部 分 在 服务 器 ( 云 ) 上 完成 ， 结 果 息 09， 娱 乐 PC， 金融 ， 社 会 科学 等 。 


灶 
吗 


通过 Web 接口 交付 给 用 户 。 众 所 周知 ， 云 计算 包含 三 种 服务 模 计算 金融 是 一 个 跨 学 科 的 领域 ， 将 数理 统计 建 模 、 分 析 和 
式 : 基础 架构 即 服务 〈IaaS )， 平 台 即 服务 (PaaS ) 和 软件 即 服 。 计算 应 用 于 定价 、 交 易 和 对 冲 决策 ， 并 分 析 这 些 决策 的 风险 。 
务 (SaaS )。 当 前 的 公开 PaaS 包括 Microsoft Windows Azure， 这 个 领域 典型 的 数理 统计 技术 包括 随机 分 析 ， 蒙 特 卡 洛 模拟 ， 


Amazon Web Services, Google App Engine 和 Sun Grid Compnute 微分 方程 ， 多 元 分 析 ， 时 间 序 列 等 。 很 多 通用 软件 兼容 了 计算 
Utility 以 及 Open Cirrus 器 等 新 兴 实 验 平 台 。Paasg 的 一 个 例子 金融 ， 例 如 MATLAB ，Mathematica，Maple 和 有 R， 另 外 还 有 专 
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门 用 于 计算 金融 的 软件 ， 例 如 Quantifi，StreamBase，Volera 和 率 的 波动 率 , r 是 利率 。 方 程 (1) 可 以 进一步 简化 为 更 简单 的 
SciFinance。 形式 ， 如 热传导 方程 。 美 式 期 权 的 左边 界 条 件 为 
Ww =K, (2.1) 
Windows Azure Node 吹 氏 苦 则 | 头 
二 | Amexote | 而 欧 氏 期 权 则 为 
HPC 一 WW=Ke ; (2.2) 
其 中 天 是 常数 ， 看 涨 期 权 的 左边 界 条 件 为 
Micrbsoft's DataCenter (IaaS/ ) | WW =0, (3.1) 
图 1 Windows Azure HPC 的 逻辑 结构 | 
Vt =0, G2) 
除 台式 电脑 外 ， 并 行 计算 还 被 广泛 应 用 于 计算 金融 领域 ， 期 权 价 格 预 测 的 最 终 条 件 为 
例如 集群 和 GPU 等 平台 。 云 并 行 计算 是 云 上 的 并 行 计算 模式 ， Vi = max(Kk — 5,0), (4) 
如 Windows Azure 的 HPC 调度 程序 ， 它 为 金融 公司 和 专业 人  ” 其 中 7 是 时 间 的 步 长 总 数 。Black-Scholes 方程 的 关键 思想 是 
员 提 ee 方 可 以 用 方程 式 中 的 “正确 ”价格 来 对 冲 资产 以 避免 所 有 的 风险 ， 
期 权 是 一 种 能 在 未 来 某 特 定时 间 以 特定 价格 买 入 或 卖 出 一 。 ”所 以 定价 的 计算 速度 很 重要 。Black-Scholes 方程 的 求解 有 两 种 
定数 量 的 某 种 特定 商品 的 权利 。 期 权 可 以 通过 Black-Scholes 方 。 方法 : 分 析 方 法 和 数值 方法 。 文献 [21] 中 ,Martin 等 人 开发 了 近 
程 或 随机 波动 模型 进行 定价 ， 并 通过 二 项 式 树 、 蒙 特 卡 洛 模拟 。” 似 解析 解 ， 而 Black-Scholes 方程 的 数值 解法 包括 有 限 差分 法 
和 有 限 差分 法 等 数值 方法 求解 。 期 权 定价 的 计算 速度 是 交易 双 ”2 下 、 蒙 特 卡 洛 模拟 B4291、 路 径 积 分 P73、 机 器 学 习 R 引 和 等。 本文 
方 取得 成 功 的 关键 因素 之 一 ， 云 并 行 计算 为 满足 期 权 定价 的 速 ” 用 有 限 差分 法 求解 Black-Scholes 方程 。 
度 需 求 提 供 了 新 的 平台 。 本 文 在 Windows Azure 的 基础 上 ， 利 
用 Black-Scholes 方程 模拟 期 权 定价 , 并 且 开 发 了 基于 有 限 差分 
软件 的 云 并 行 计算 AzureOP 进行 期 权 定价 。 


1 ”AzureOP 方法 


首先 介绍 如 何在 云 上 使 用 AzureOP 进行 期 权 定价 。 然 后 介 
绍 如 何 设计 AzureOP: 包括 三 种 Black-Scholes 方程 的 数值 求解 
方案 ， 这 三 种 方案 的 云 并 行 实现 ， 以 及 云 上 使 用 到 的 软件 
且 在 第 一 次 运行 之 前 ， 将 软件 部 署 到 Windows Azure。 最 后 介 
绍 AzureOP 的 计算 结果 。 
1.1 软件 
于 AzureOP 是 基于 云 平台 Azure 的 软件 ， 安 装 过 程 由 平 
台 配 置 过 程 组 成 ， 有 具体 细节 将 在 后 面 讨 论 。 在 成 功 部 署 之 后 ， 
除了 超级 计算 和 背后 的 大 数据 存储 的 功能 ，AzureOP 的 使 用 类 
似 于 任何 桌面 软件 。 
AzureOP 的 GUI 如 图 2 所 示 ，AzureOP 的 界面 由 三 部 分 组 
成 :“ 选 择 选项 *, “参数 设置 "和 命令 按钮 。 选 项 样式 可 以 从 
“option selection” 部 分 中 选择 ， 参 数 可 以 从 “parameter setting”* 冲 
分 进行 设置 。 之 后 ， 按 下 “start” 按 钮 启动 AzureOP。 AzureOP 
的 中 间 结 果 存 储 在 Azure blob 中 ， 按 下 ‘下 etch Results” 按 钮 可 以 
将 计算 结果 提取 到 AzureOP 的 输出 帧 。 


1.2 ”模型 

AzureOP 通过 求解 Black-Scholes 方程 来 做 期 权 定价 。 
Black-Scholes 模型 是 一 个 二 阶 偏 微分 方程 ， 它 描述 了 期 权 定价 
从 终止 条 件 递归 到 价格 为 零 的 结束 时 期 的 过 程 。 从 原始 文献 1 


推导 出 ee 方程 后 ， 得 到 : 
+162522Y +rST ry =0, (1) 


轴 权 价格 ， S 是 股票 价格 ，t 是 时 间 ，5 是 股票 收益 


其 中 : 是 
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1.3 显 式 方 案 


AzureOP 由 三 个 方案 组 成 ， 


2 ”AzureOP 的 界面 


程 : 显 式 方案 、 隐 式 方案 和 Crank-Nicolson 方案 。 


案 离 散 Black-Scholes 方程 ， 本 文 定 义 下 面 


ee 


OV 加 Ti+1L 一 


ot 


az _ 
0S2 


aV 
GAY 


» 


At 


1 1 1 
Va -2 + 


(AS)? 


_ Vtl ynt! 
2AS 


用 于 数值 求解 Black-Scholes 方 


过 显 式 方 


(5.1) 


(5.2) 


(5.3) 


其 中 ; n 是 时 间 步 长 的 计数 器 , i 是 网 格 (股票 价格 ) 的 计数 器 。 
将 这 些 算 子 (5.1-5.3) 代入 方程 (1) 得 到 
Vr = AiVit! + BiVit! 十 CT 所 


其 中 ， 


Bi = 


Su | ee 
i 1+rAt 


1 
1 一 
1+rAt ( 


2 ) 
2AS 己 
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(As)2 己 
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(6.3) 
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1 一 二 GG | (64) 
方程 6.1) 可 以 得 到 ， 期 权 价格 从 最 终 状 态 (4) 向 当前 
时 间 步 长 二 后 向 递归 演变 ， 每 个 网 格 〈 股 票 价格 ) 的 期 权 价格 
从 1 到 工 1 进行 逐一 计算 ， 其 中 了 是 网 格 步 长 的 总 数 。 期 权 价 
格 预测 的 左边 界 条 件 可 以 被 离散 化 为 : 
V8 = Ke nat), (7.1a) 
或 者 
Ww =0, (7.1b) 
期 权 价 格 方程 (3.2〉 的 右边 界 条 件 可 以 被 离散 化 为 : 
Vr = (7.2) 


其 中 6 和 不 需要 计算 ， 可 以 直接 用 边界 条 件 代 栓 。 另 外 ， 为 
了 保持 显 式 方案 的 稳定 性 ， 应 该 满足 条 件 At < 和。 


1.4 隐 式 方案 
AzureOP 还 可 以 通过 隐 式 方案 离散 Black-Scholes 方程 , 算 
子 定义 如 下 : 
OV Vi-2W + 
O52 (AS)?  ， (8.1) 
oav Vi 一 Vi 
- -- 一 (2 


将 方程 (5.1), (8.1) 和 (8.2) 代 入 Black-Scholes 方程 (1)， 
可 以 得 到 : 


Vi = AiVi + BiVr + CV (9.1) 
其 中 : 
At 6252 srs 
和 = = (AS)? 本 人 02) 
B=1+55 4rAt 93 
(AS)2 rat, ( * ) 
B252 rs 
4 = = 二 本 | 94 
将 方程 (9.1) 转换 成 矩阵 形式 ， 可 以 得 到 
BC 0 入 0 0 0 ww Vit1 — AiVe hh 
4 B2 C * 0 0 0 Va 慑 +1 E 
;  ” : 1: : |= : =| : |.(10) 
0 0 0 4-z Bz Czllviz vi! ti-2 
000 “= 0 4a Bl Lyre wn) be 


其 中 记 是 线性 系统 方程 (10) 的 等 式 右 侧 部 分 。 为 了 求解 系统 
方程 (10),， 使 用 三 对 角 和 矩阵 算法 (TDMA)。 与 显 式 方案 相 比 ， 
隐 式 方案 在 任何 At 和 As 情况 下 都 更 稳定 。 

另外 , 由 于 线性 系统 方程 (10) 的 系数 矩阵 有 三 对 角 线 ( 左 
对 角 线 , 中 对 角 线 和 右 对 角 线 ) 所 以 这 些 对 角 线 的 信息 被 压缩 
成 三 个 矢量 〈 左 矢量 ， 中 矢量 和 右 矢 量 )， 并 且 忽 略 所 有 的 零 。 
因此 ， 在 一 个 7z 的 系数 矩阵 中 ， 左 矢量 的 大 小 为 天 1， 中 矢量 
的 大 小 为 也 右 矢 量 的 大 小 为 大 1。 
1.5 ”Crank-Nicolson 方案 

AzureOP 求解 Black-Scholes 方程 的 第 三 种 内 置 方案 是 
Crank-Nicolson 方案 ， 算 子 定 义 如 下 : 


ov — (V2 H+ -2 十 由 
952 2(AS)? 


(11.1) 
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-人 二 一 + ) (11.2) 
os 4AS 
n+l1 n 
V= 0 (11.3) 


将 方程 (5.1) (11.1) (11.2) 和 (11.3) 代入 Black-Scholes 
方程 (1)， 可 以 得 到 : 


AiVic1 + BiVi + CiVits = fo (12.1) 
其 中 : 
Ai = (262i2 —iri) (12.2) 
B, = -262 2r) (12.3) 
c=(: 62i2 + iri) (12.4) 


f=(-i62i ve +(- 志 +362i2 ++37) Vt1 十 
(~#6222 — tri) Ved (12.5) 
将 方程 (12.1) 转换 成 矩阵 形式 ， 可 以 得 到 : 


BI CGC 0 … 0 0 0 WA f+ — AiVe 局 

4 B, G 汪 0 0 0 V2 jt b 

000 hs Bos cm Hi “|,, 

0 0 0 0 Ar1 Bi-JLri frtl— Cv Fi 

(13) 

其 中 Fi 是 线性 系统 方程 (13 ) 的 等 式 右 侧 部 分 , 并 且 采 用 TDMA 
来 求解 该 系统 。 
1.6 云 并 行 计算 实现 


AzureOP 由 上 述 三 种 有 限 差分 方案 的 云 并 行 计算 实现 组 成 ， 
使 用 的 集成 开发 环境 是 MS Visual Studio 2011 (MS-VS)。 

为 了 在 尽 可 能 低 的 成 本 下 实现 快速 收敛 ,使 用 MS-MPTI 库 
对 Black-Scholes 方程 的 三 种 有 限 差分 方案 进行 并 行 化 处 理 ， 
在 微软 的 Windows Azure 上 运行 ， 以 获得 可 扩展 的 实例 数量 及 
满足 用 户 对 速度 的 需求 。 
在 MS-VS 上 开发 并 行 代 码 ， 首 先 需 要 安装 MPI 库 ， 本 文 
选择 的 MPI 库 是 HPC Pack 2008 R2 MS-MPI Redistributable 
Package。 安 装 完成 后 ， 重 复 检 查 MS-VS 以 确保 HPCPACK 库 
被 正确 引用 。 然 后 , 利用 并 行 库 MS-VS 开发 了 
行 隐 式 方案 和 并 行 Crank-Nicolson 方案 。 

在 显 式 方 案 的 方程 6.1) 中， 求解 当前 时 间 步 长 Vr 的 期 
权 价 格 需要 后 续 三 个 时 间 步 长 的 值 : 严 国 2， 了 友和 Fr。 在 
求解 当前 时 间 的 期 权 价格 时 ， 由 于 后 续 时 间 步 长 的 期 权 价格 的 
所 有 值 都 是 可 用 的 ， 所 以 这 些 值 可 以 同时 求解 。 因 此 本 文 设计 
了 一 个 如 图 3 和 4 所 示 的 云 并 行 显 式 方案 。 

图 3 描述 了 显 式 方 案 的 云 并 行 算 法 。 如 图 3 所 示 ， 每 个 时 
间 步 长 中 期 权 价 格 的 值 是 同时 求解 的 ， 而 不 是 逐一 求解 ， 并 且 
是 在 当前 时 间 步 长 到 下 一 个 时 间 步 长 之 间 的 时 间 间 隔 中 进行 求 
解 。 
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行 显 式 方案 、 


4 描述 了 AzureOP 显 式 方案 中 的 实现 。 如 图 4 所 示 ， 
MS-MPI 进程 分 为 三 类 : 主 服务 器 《初始 进程 )， 网 格 服务 器 和 
工作 节点 。 主 服务 器 〈 初 始 进程 ) 负责 存储 和 更 新 期 权 定价 的 
所 有 值 ， 网 格 服务 器 负责 计算 每 个 网 格 的 期 权 定价 ， 节 点 工作 
人 员 则 负责 计算 方程 6.2) 中 的 省 ， 方 程 〈6.3) 中 的 B; 和 方 


(nN | nn 
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程 (6.4) 中 的 Ci。 


对 


5 隐 式 方案 的 云 并 行 算法 《算法 2) 


图 3 显 式 方案 的 云 并 行 算法 (算法 1) 


The Grid The Grid 
ASPNET/ - ServerforF; ServerforF; 
CGI User Price 
到 Interface The Grid 
台 Server 日 
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forGof 
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Web Role Storage 


Windows Azure Worker Role 


6 AzureOP 隐 式 方案 的 实现 


图 4 ” AzureOP 显 式 方 案 的 实现 


图 6 描述 了 AzureOP 隐 式 方案 的 实现 。 从 图 6 可 以 看 到 ， 


为 了 更 加 清晰 地 表达 ， 图 3 和 4 中 讨论 的 云 并 行 显 式 方案 ” MS-MPI 进程 分 为 三 类 : 主 服 务 器 〈 发 起 者 )， 网 格 服务 器 和 工 
也 可 以 如 算法 1 所 述 : 作 节 点 。 主 服务 器 负责 存储 和 更 新 期 权 价格 的 值 , 组 装 方程 (10) 
算法 1 基于 云 并 行 计算 的 显 式 方案 。 的 系数 矩阵 和 其 右 侧 矢量 F;， 并 求解 线性 系统 。 网 格 服务 器 负 
a) 发 起 方 〈 主 服务 器 ) 将 终止 条 件 (4) 初始 化 为 Vin+1; 责 计算 方程 10) 的 右 侧 矢量 F， 并 收集 方程 (9.2) 的 4i, 方 
b) 主 服务 器 将 方程 (7.1) 和 (7.2) 的 边界 条 件 设置 为 Vin+3; 程 (9.3) 的 Bi 或 方程 (9.4) 的 Ci。 节 点 工作 人 员 则 负责 计算 
二 c) 在 每 一 个 时 间 步 长 中 : 方程 (9.2) 的 4;,， 方程 (9.3) 的 Bi 或 方程 (9.4) 的 Ci 的 值 。 
(a) 网 格 服务 器 初始 化 Vi"; 云 并 行 计算 的 隐 式 方案 也 可 以 如 算法 2 所 述 : 
(b) 计算 方程 (6.2) - (6.4) 中 的 Ai，Bi 和 Ci， 随 后 工作 节点 将 这 些 算法 2 基于 云 并 行 计算 的 隐 式 方案 。 
值 发 送 给 对 应 的 网 格 服务 器 ; a) 发 起 方 〈 主 服务 器 ) 将 终止 条 件 〈4) 初始 化 为 Vin+1; 
(c) 网 格 服务 器 计算 Vi"， 并 发 送 给 主 服务 器 ; b) 主 服务 器 将 方程 (7.1) 和 “(7.2) 的 边界 条 件 设 置 为 Vin+1; 
(Cd) 主 服 务 器 将 vi" 保存 到 Windows Azure Storage 上 的 价格 容器 ; c) 在 每 一 个 时 间 步 长 中 : 
(e) 主 服务 器 根据 Vi" 更 新 Vi"*1; (a) 网 格 服务 器 初始 化 Vi"; 
结束 (b) 计算 方程 (9.2) ~ (9.4) 中 的 Ai，Bi 和 Ci， 随后 工作 节点 将 这 些 
在 开发 了 算法 1 的 云 并 行 显 式 方案 后 ， 本 文 开发 了 云 并 行 值 发 送 给 对 应 的 网 格 服务 器 ; 
殷 式 方案 。 在 隐 式 方案 的 线性 系统 方程 (10〉 中， 通过 求解 系 (c) 网 格 服务 器 计算 方程 (16) 中 的 Fi 值 ， 并 发 送 给 主 服 务 器 ; 
统 可 以 同时 得 到 期 权 价 格 Vr 的 所 有 值 ,为 了 加 速 这 一 过 程 , 在 (d) 主 服务 器 组 装 系数 矩阵 和 右 侧 矢量 Fi, 并 求解 线性 系统 方程 (10); 
显 式 方案 中 开发 的 云 并 行 实现 策略 的 基础 上 ， 可 以 计算 出 工作 (e) 主 服务 器 将 Vi" 保 存 到 Windows Azure Storage 上 的 价格 容器 ; 
节点 中 的 系数 和 矩阵， 并 解决 主 服务 器 中 的 线性 系统 问题 。 因 此 (f) 主 服务 器 根据 vi" 更 新 Vin*3; 
本 文 设计 了 一 个 如 图 5 和 6 所 示 的 云 并 行 隐 式 方案 。 结束 
图 5 描述 了 隐 式 方案 的 云 并 行 算法 。 如 图 5 所 示 ， 利 用 图 云 并 行 算法 的 Crank-Nicolson 方案 也 可 以 类 比 图 5 的 描述 。 


6 中 列 出 的 服务 器 ， 可 以 并 行 地 计算 和 组 装 系数 矩阵 ， 当 前 时 ”如 图 5 所 示 ， 与 隐 式 方案 类 似 ， 方程 (13)〉 中 的 系数 矩阵 被 同 
间 步 长 的 期 权 价格 的 所 有 值 也 可 以 由 线性 系统 方程 (10) 求 解 ， 时 计算 和 组 装 ， 当 前 时 间 步 长 的 所 有 期 权 价 格 的 值 由 线性 系统 
并 且 此 过 程 一 直 持续 。 方程 (13) 进行 求解 ， 并 且 这 个 计算 过 程 在 所 有 时 间 步 长 中 持 
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续 步骤 
图 7 描述 了 AzureOP 的 Crank-Nicolson 方案 的 实现 。 主 服 ee 

ds i 2 ”计算 结果 
务 器 负责 存储 和 更 新 期 权 价格 的 值 ,组 装 系数 矩阵 和 右 侧 矢量 ， 
并 求解 线性 系统 方程 (13)。 网 格 服 务 器 负责 计算 右 侧 矢量 的 Fi 本 章 将 介绍 使 用 AzureOP 进行 期 权 定 价 的 计算 结果 ,并 将 
值 ， 并 收集 方程 (12.2) 的 4i, 方程 (12.3) 的 Bi, 方程 (12.4) 显 式 方案 ， 隐 式 方案 和 Crank-Nicolson 方案 的 结果 一 一 列 出 。 
的 Ci 或 方程 (12.5) 的 的 值 。 如 图 7 所 示 ， 与 图 6 中 实现 的 在 显 式 方案 中 , 总 时 间 为 0.01, 离散 化 为 230、500 或 1000 
总 式 方案 相同 , 图 7 中 的 工作 节点 负责 计算 方程 (12.2) 的 4， 步 ， 股 价 St 为 1， 离散 化 为 23、50 或 100。AzureOP 的 时 间 成 


方程 (12.3) 的 Bi, 方程 (12.4) 的 Ci 或 方程 (12.5) 的 的 值 。 本 和 费用 如 表 1 所 示 ; 时 间 步 长 = 250、 网 格 步 长 = 100 (7 = 
100) 时 的 定价 曲面 计算 结果 如 图 8 所 示 ; 显 式 方案 在 t=1.00， 

时 间 步 长 = 250 和 7 了 = 100 时 的 定价 曲线 如 图 9 所 示 ， 显 式 方 

案 在 s=0.49 〈s 是 股票 价格 的 当前 值 )， 时 间 步 长 =250 和 7= 

100 时 的 定价 曲线 如 图 10 所 示 。 

| Sen 在 隐 式 方案 中 ， 总 时 间 为 1， 离 散 化 为 250、500 或 1000 

步 ， 股 票 价格 St 为 1， 离散 化 为 25、50 或 100。AzureOP 的 时 


间 成 本 和 费用 如 表 2 所 示 ; 时 间 步 长 = 250、 网 格 步 长 = 100 
ee es a 3 四 A ZE 本 a 
et oe oe oy (T= 100) 时 的 定价 曲面 计算 结果 如 图 11 所 示 ; 隐 式 方案 在 1 
Windows Azure Worker Role = 1.00,， 时 间 步 长 =250 和 1= 100 时 的 定价 线 如 图 12 所 示 ; 
隐 式 方案 在 s = 0.08， 时 间 步 长 = 250 和 7Z= 100 时 的 定价 曲线 


图 7 AzureOP 的 Crank-Nicolson 方案 的 实现 


如 图 13 所 示 。 


云 并 行 算法 的 Crank-Nicolson 方案 也 可 以 如 算法 3 所 述 : 在 Crank-Nicolson 方案 中 , 总 时 间 为 1, 离散 化 为 250、500 
算法 3 基于 云 并 行 计算 的 Crank-Nicolson 方案 。 或 1000 步 ， 股 价 St 为 1， 离 散 化 为 25、50 或 100。AzureOP 
a) 发 起 方 ( 主 服务 器 ) 将 终止 条 件 (4) 初始 化 为 Vin+1; 的 时 间 成 本 和 费用 如 表 3 所 示 ; 时 间 步 长 = 250、 网 格 步 长 = 
b) 主 服务 器 将 方程 (7.1) 和 (7.2) 的 边界 条 件 设置 为 Vin+1; 100 (7 = 100) 时 的 定价 曲面 计算 结果 如 图 14 所 示 ; Crank- 
c) 在 每 一 个 时 间 步 长 中 : Nicolson 方案 在 += 1.00， 时 间 步 长 = 250 和 7= 100 时 的 定价 

(a) 主 服务 器 初始 化 win; 线 如 图 15 所 示 ; Crank-Nicolson 方案 在 s = 0.17， 时 间 步 长 
(b) 计算 方程 (12.2) - (12.5) 中 的 Ai，Bi，Ci 和 天 ,随后 工作 节点 。 =250 和 7= 100 时 的 定价 曲线 如 图 16 所 示 。 

将 这 些 值 发 送 给 对 应 的 网 格 服务 器 ; 本 文选 择 的 计算 实例 规模 非常 小 ， 一 共 购 买 了 25 个 实例 ， 
(c) 网 格 服务 器 计算 方程 (13) 中 的 户 值 ， 并 发 送 给 主 服务 器 ; 且 共 享 CPU 内 核 和 768 MB 内 存 ， 每 小 时 收费 0.02 美元 。 
《d) 主 服务 器 组 装 系数 矩阵 和 右 侧 矢量 fi, 并 求解 线性 系统 方程 (13); 表 1 是 AzureOP 选择 显 式 方案 的 计算 结果 。 从 表 1 可 以 得 
Ce) 主 服务 器 将 Vi" 保 存 到 Windows Azure Storage 上 的 价格 容器 ; 到 ， 显 式 方案 的 时 间 成 本 从 几 分 钟 到 几 小 时 不 等 ， 而 计算 的 开 
(f) 主 服 务 器 根据 Vi" 更 新 Vin'3; 销 很 低 。 因 此 在 保证 最 低 费 用 的 情况 下 ， 可 以 选择 一 个 最 佳 的 
结束 参数 时 间 步 长 和 网 格 步 长 ， 以 获得 良好 的 模拟 结果 ， 表 1 的 最 

1.7 Windows Azure 的 安装 多 选 择 是 时 间 步 长 =1000、7= 50。 
成 功 安 装 AzureOP， 首 先 需要 一 个 Windows Azure 帐户 。 图 8 是 Azure 显 式 方案 的 计算 结果 。 从 图 8 可 以 看 到 ， 随 
之 后 ,使 用 Windows Azure HPC 调度 程序 来 提交 、 管 理 和 获取 着 时 间 步 长 的 增加 ， 期 权 价格 几乎 保持 不 变 ， 出 现 这 种 情况 的 


三 种 方案 的 代码 结果 ， 其 中 调度 程序 包括 工作 节点 〈 头 节点 ， 原因 是 ， 为 了 满足 显 式 方案 稳定 性 的 要 求 ， 本 文选 择 了 一 个 非 
计算 节点 和 前 端 节 点 ) 和 数据 存储 (Windows Azure Storage)。 ” 常 小 的 At， 约 等 于 105， 这 导致 了 一 个 非常 小 的 时 间 间 隔 了 = 


Windows Azure HPC 调度 程序 将 Windows Azure 转变 为 基于 云 ”0.01， 因 此 显 式 方案 几乎 没有 进展 。 

的 超级 计算 机 。 所 有 步骤 都 需要 通过 在 AppConfigure 应 用 程序 表 1 AzureOP 显 式 方案 的 时 间 花 销 和 费 

中 进行 配置 并 填写 Windows Azure 帐户 、 管 理 员 帐户 以 及 头 节 时 间 网 格 步 长 -25 网 格 步 长 =50 网 格 步 长 =100 

点 、 计 算 节 点 、 前 端 节点 的 数量 等 信息 来 完成 的 。 运 行 命令 步 长 时 间 花 销 /s 费用 /美元 时 间 花 销 费用 /美元 “时间 花 销 A。 费用 /美元 

hpcfiwutil 可 以 打开 计算 机 节点 间 通 信 的 防火 墙 配置 。 250 55.60 0.10 202.44 0.10 1096.21 0.10 
Windows Azure HPC 调度 程序 准备 就 绪 后 ， 远 程 桌面 连接 500 219.37 0.10 791.70 0.10 4339.11 0.20 

即 可 访问 计算 节点 实例 。 之 后 ， 开 发 的 代码 可 以 通过 Windows 1000 877.07 0.10 3145.75 0.10 17319.28 0.50 

Azure Web Role 的 图 形 用 户 界面 提交 给 Windows Azure。 第 一 


次 安装 之 后 , 后续 AzureOP 的 使 用 都 不 需要 涉及 本 节 中 的 所 有 
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图 8 时间 步 长 =250， 网 格 步 长 =100 (7=0.01, %=1.00, 6=0.30, 一 0.01， 
K=0.50) 时 ， 显 式 方 案 的 定价 曲面 


为 了 比较 不 同时 间 步 长 的 计算 结果 ， 本 文 绘制 了 t= 1.00， 
以 及 时 间 步 长 分 别 为 250、500 和 1000 时 的 期 权 价格 数值 结果 ， 
如 图 9 所 示 。 
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图 9 ”时 间 步 长 =250， 网 格 步 长 =100 (f=1.00，7=0.01, S=1.00, 6=0.30， 
三 0.01, K=0.50) 时 ， 显 式 方案 的 定价 曲线 


图 9 哲 


绘 了 显 式 方案 的 期 权 价 格 数值 结果 。 从 图 9 中 可 以 
看 到 ， 和 迭代 次 数 等 于 250、500 和 1000 的 期 权 价格 数值 结果 完 
全 一 致 。 时 间 步 长 的 选择 不 会 影响 期 权 价格 预测 数值 结果 的 精 
度 ， 但 会 影响 数值 结果 的 解析 度 。 

为 了 研究 期 权 价格 随时 间 进 度 产 生 的 变化 ， 本 文 绘制 了 s 
= 0.49， 以 及 时 间 步 长 分 别 为 50、500 和 1000 时 的 期 权 价格 
数值 结果 ， 如 图 10 所 示 。 
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图 10 时 间 步 长 =250, 股 票 当 前 价格 =0.49, 网 格 步 长 =100(T=0.01, 
St=1.00, 6=0.30, r=0.01，K=0.50) 时 , 显 式 方案 的 定价 曲线 


从 图 10 可 以 看 出 , 迭代 次 数 等 于 250、500 和 1000 的 期 权 
价格 数值 结果 完全 一 致 ， 说 明 时 间 步 长 的 选择 不 影响 期 权 价格 
的 数值 结果 。 然 而 ， 随 着 时 间 的 推移 ， 期 权 价格 几乎 没有 发 生 
变化 ， 原 因 正 如 上 文 讨论 的 。 


t 


案 在 s = 0.08 时 的 期 权 价格 
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在 分 析 显 式 方案 的 结果 之 后 ， 看 一 下 隐 式 方案 的 结果 ， 如 
表 2、 图 11~13 所 示 。 表 2 是 隐 式 方案 的 计算 结果 。 从 表 2 可 
以 看 出 ， 虽 然 隐 式 方案 的 时 间 成 本 从 几 分 钟 到 几 小 时 不 等 ， 但 
时 间 步 长 和 网 格 步 长 每 个 组 合 的 开销 都 很 低 。 参 数 的 最 佳 选择 
是 时 间 步 长 = 1000，7= 50 和 时 间 步 长 =500，7= 100。 


表 2 AzureOP 隐 式 方案 的 时 间 花 销 和 费用 
时 间 网 格 步 长 =25 网 格 步 长 =50 网 格 步 长 =100 
步 长 ”时间 花 销 /s ”费用 /美元 ”时 间 花 销 /s ”费用 /美元 ”时间 花 销 /s ”费用 /美元 
250 50.21 0.10 168.20 0.10 655.56 0.10 
500 199.07 0.10 660.41 0.10 2600.02 0.10 
1000 786.28 0.10 2607.49 0.10 10282.66 0.30 


比较 表 1 和 2 可 以 得 到 ， 几 乎 每 个 时 间 步 长 和 网 格 步 长 的 
组 合 中 ， 隐 式 方 案 都 比 显 式 方案 快 。 当 期 权 到 期 时 ， 期 权 价格 
无 限 趋 近 于 0。 

图 11 是 隐 式 方案 的 计算 结果 。 从 图 11 可 以 看 出 ，7= 0.01 
的 时 间 间 隔 足 够 长 ， 因 此 随 着 时 间 的 推进 ， 当 期 权 到 期 时 ， 期 
权 价 格 接近 零 。 比 较 图 8 和 11 可 以 得 到 ,图 11 比 图 8 更 平滑 ， 
原因 在 于 隐 式 方案 选 定 的 时 间 间 隔 了 远大 于 显 式 方 案 的 时 间 间 


隔 。 


图 11 


时 间 步 长 =250， 网 格 步 长 =100 (7=0.01, $=1.00, 6=0.30, :=0.01, 
K=0.50) 时 ， 隐 式 方案 的 定价 曲面 


图 12 绘制 了 隐 式 方案 的 期 权 价格 数值 结果 。 从 图 12 可 以 
看 到 ， 和 迭代 次 数 等 于 250、500 和 1000 的 期 权 价格 数值 结果 完 
全 一 致 。 时 间 步 长 的 选择 不 会 影响 期 权 价格 预测 数值 结果 的 精 
度 ， 但 会 影响 数值 结果 的 解析 度 。 
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图 12 时间 步 长 =230， 网 格 步 长 =100 (于 1.00, 7=1.00, S=1.00, 6=0.30， 
/0.01, K=0.50) 时 ， 隐 式 方案 的 定价 曲线 


为 了 研究 随时 间 变 化 的 特定 股票 价格 ， 本 文 绘制 了 隐 式 方 
线 ， 如 图 13 所 示 。 从 图 13 可 以 
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看 出 ， 和 迭代 次 数 等 于 250、500 和 1000 的 期 权 价格 数值 结果 完 
全 一 致 。 与 显 式 方 案 相 同 ， 时 间 步 长 的 选择 不 会 影响 期 权 价 格 
预测 数值 结果 的 精度 ， 但 会 影响 数值 结果 的 解析 度 。 


Option Price 


10 


图 13 


时 间 步 长 =250， 股 票 当 前 价格 =0.08， 网 格 步 长 =100 (7=10.00， 
SF 六 1.00, 6=0.30, /=0.01, K=0.50) 时 ， 隐 式 方案 的 定价 曲线 


从 图 13 还 可 以 看 出 , 迭代 次 数 等 于 250、500 和 1000 的 期 
权 价格 数值 结果 完全 一 致 ， 这 意味 着 时 间 步 长 的 选择 不 会 影响 
期 权 价格 预测 数值 结果 的 精度 ， 但 会 影响 数值 结果 的 解析 度 。 

最 后 ， 将 计算 结果 绘制 出 来 ， 如 表 3、 图 14~16 所 示 。 表 
3 是 Crank-Nicolson 方案 的 计算 结果 。 从 表 3 可 以 看 到 ，Crankc- 
Nicolson 方案 的 时 间 成 本 从 几 分 钟 到 几 小 时 不 等 , 费用 也 很 低 。 
时 间 步 长 和 网 格 步 长 参数 的 最 佳 组 合 是 时 间 步 长 = 1000， 网 格 
步 长 = 50 和 时 间 步 长 = 500， 网 格 步 长 = 100。 


表 3 AzureOP 的 Crank-Nicolson 方案 的 时 间 花 销 和 费用 


时 间 网 格 步 长 =25 网 格 步 长 =50 网 格 步 长 =100 


步 长 ”时间 花 销 /s ”费用 /美元 ”时间 花 销 /s ”费用 /美元 ”时 间 花 销 /s ”费用 /美元 


250 57.87 0.10 198.27 0.10 782.11 0.10 
500 228.35 0.10 788.86 0.10 3112.62 0.10 
1000 905.65 0.10 3125.20 0.10 12418.10 0.30 


图 14 是 具有 特定 参数 的 Crank-Nicolson 方案 的 定价 曲面 ， 
从 图 14 可 以 看 到 ， 随 着 时 间 的 推移 ， 当 期 权 到 期 时 ， 期 权 价格 
无 限 趋 近 于 零 。 
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图 14 时 间 步 长 =250， 网 格 步 长 =100 (7=1.00, % =1.00, 6=0.30， 
/=0.01, K=0.50) 时 ，Crank-Nicolson 方案 的 定价 曲面 


图 15 绘制 了 1= 1.00 时 Crank-Nicolson 方案 的 期 权 价格 曲 


格 曲线 。 从 图 16 可 以 看 到 ,和 迭代 次 数 等 于 250、500 和 1000 的 
期 权 价格 数值 结果 完全 一 致 ， 这 意味 着 价格 曲线 对 时 间 步 长 
不 敏感 。 从 图 16 还 可 以 看 到 ， 随 着 时 间 的 推进 ， 期 权 价 格 趋 近 
于 零 ， 这 与 先前 的 假设 相 吻 合 。 


Stock Price 


时 间 步 长 =250， 网 格 步 长 =100 (f=1.00, 7=1.00, S=1.00, 6=0.30， 


图 15 


/一 0.01, K=0.50) 时 ，Crank-Nicolson 方案 的 定价 曲线 


Tin 


图 16 ”时间 步 长 =250, 股 票 当前 价格 =0.17, 网 格 步 长 =100 ( 夺 1.00， 
7=1.00, S=1.00, 5=0.30, /=0.01, K=0.50) 时 ，Crank-Nicolson 方案 的 定 


价 曲线 
3 ”结束 语 
本 文 利 用 Windows Azure 的 HPC 调度 程序 ， 开 发 了 基于 


云 并 行 计算 的 期 权 定价 方案 ， 该 方案 以 较 低 的 费用 提供 了 低 风 
险 和 高 速度 ， 并 给 出 了 两 种 方案 的 期 权 定价 仿真 结果 。 
正如 上 述 2.3 节 中 讨论 的 那样 ， 主 服务 器 (发 起 者 ) 利用 
TDMA 方法 求解 线性 系统 方程 (10)。 为 了 进一步 加 快 求解 线 
性 系统 的 计算 速度 ，TDMA 可 以 通过 并 行 计算 来 实现 ， 而 现 
的 并 行 TDMA 研究 包括 集群 和 GPU。 文献 [29] 中 , Swarztrauber 
等 人 开发 了 一 个 通用 三 对 角 和 矩阵 的 并 行 算法 ， 文 献 [30] 中 ， 
Nathan 等 人 开发 了 一 个 求解 三 对 角 和 矩阵 的 并 行 算法 ， 文 献 [31] 
中 ，Lin 等 人 则 开发 了 一 个 基于 图 的 三 对 角 和 矩阵 算法 。 

期 权 定价 的 研究 中 ， 美 式 期 权 的 定价 问题 ， 通 常 采用 各 种 
偏 微分 方程 数值 方法 来 求解 。 针 对 此 类 计算 量 大 、 计 算 时 间 长 
的 问题 ,并行 计 算 提 供 了 有 效 的 方法 B233]。 计 算 速 度 固 然 重要 ， 
但 现 有 的 美式 期 权 定 价 研究 主要 将 重心 集中 在 数学 模型 的 改进 


线 。 从 图 15 可 以 看 到 ,迭代 次 数 等 于 250、500 和 1000 的 期 权 
价格 数值 结果 完全 一 致 ， 这 意味 着 时 间 步 长 的 选择 不 会 影响 期 
权 价格 预测 数值 结果 的 精度 ， 但 会 影响 数值 结果 的 解析 度 。 
图 16 绘制 了 当 s = 0.17 时 ，Crank-Nicolson 方案 的 期 权 价 


和 提高 运算 速度 上 FE 1， 很 少 可 以 兼顾 数据 存储 成 本 、 计 算 成 
本 和 硬件 开销 问题 。 而 云 计算 基于 互联 网 ,将 规模 化 的 计算 能 
存储 空间 、 开 发 平台 和 软件 服务 提供 给 用 户 ,用 户 只 需要 支付 
定 的 费用 ， 便 可 以 随时 随地 通过 互联 网 使 用 云 计算 提供 的 资源 
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[38,39], 
金融 领域 云 并 行 计算 的 一 个 关键 问题 是 应 该 选择 哪 种 部 署 
模式 : 公有 云 还 是 私有 云 ? 对 于 没有 兴趣 建设 和 维护 数据 中 心 
的 金融 公司 来 说 ， 公 有 云 是 一 个 合适 的 选择 ， 因 为 金融 公司 只 


需要 向 laaS 和 PaaS 提供 商 支付 账单 。 通 常 来 说 ， 公 有 云 的 成 
本 比 私有 云 低 , 因为 公共 数据 中 心 的 规模 往往 比 私有 云 大 得 多 。 
此 外 ， 由 于 公共 Paas 的 开发 人 员 往往 来 自 领先 的 IT 公司 ， 因 


此 公共 Paag 在 技术 上 也 比 私 有 云 更 先进 。Windows Azure 此 类 
云 上 的 并 行 计算 模式 ， 为 金融 公司 和 专业 人 员 提 供 了 可 扩展 和 
低 成 本 的 并 行 计算 ， 也 为 满足 美式 期 权 定价 的 速度 需求 提供 了 
新 的 平台 。 而 本 文 提供 已 成 型 的 试点 云 软件 AzureOP， 将 计算 
和 存储 都 进行 了 交付 ,大 大 提高 了 计算 效率 ,缩短 了 执行 时 间 ， 

相对 于 自己 购买 各 种 资源 来 说 ， 以 较 低 的 成 本 完成 了 复杂 的 、 

大 量 的 计算 任务 ， 简 单 高 效 ， 且 使 用 方法 简便 快捷 。 私 有 云 则 
由 公司 的 敏感 数据 提供 了 更 多 的 安全 保障 ， 未 来 将 在 这 一 
领域 开展 更 多 的 项 目 。 
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